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Clumpy structure in the debris disk around Vega has been previously reported at millimeter wave- 
lengths and attributed to concentrations of dust grains trapped in resonances with an unseen planet. 
However, recent imaging at similar wavelengths with higher sensitivity has disputed the observed 
structure. We present three new millimeter-wavelength observations that help to resolve the puzzling 
and contradictory observations. We have observed the Vega system with the Submillimeter Array 
(SMA) at a wavelength of 880 /im and angular resolution of 5" ; with the Combined Array for Re- 
search in Millimeter- wave Astronomy (CARMA) at a wavelength of 1.3 mm and angular resolution 
of 5"; and with the Green Bank Telescope (GBT) at a wavelength of 3.3 mm and angular resolution 
of 10" . Despite high sensitivity and short baselines, we do not detect the Vega debris disk in either 
of the intcrferomctric data sets (SMA and CARMA), which should be sensitive at high significance 
to clumpy structure based on previously reported observations. We obtain a marginal (3cr) detection 
of disk emission in the GBT data; the spatial distribution of the emission is not well constrained. 
We analyze the observations in the context of several different models, demonstrating that the ob- 
servations are consistent with a smooth, broad, axisymmetric disk with inner radius 20-100 AU and 
width > 50 AU. The intcrferomctric data require that at least half of the 860 ^m emission detected by 
previous single-dish observations with the James Clerk Maxwell Telescope be distributed axisymmet- 
rically, ruling out strong contributions from flux concentrations on spatial scales of <100AU. These 
observations support recent results from the Plateau de Bure Interferometer indicating that previous 
detections of clumpy structure in the Vega debris disk were spurious. 

Subject headings: circumstcllar matter — planetary systems — planet-disk interactions — stars: in- 
dividual (Vega) 



1. INTRODUCTION 

The presence of tenuous, second generation dust disks 
around main sequence stars came as a surprise after the 
launch of the Infrared Astronomical Satellite (IRAS) in 
1983 (Aumann et al. 1984). A variety of physical pro- 
cesses is expected to remove orbiting dust grains on 
timescales shorter than the stellar age, primarily through 
collisional grinding and subsequent ejection of small 
grains by stellar radiation (e.g., Wyatt 2005; Strubbe & 
Chiang 2006). It is therefore thought that the circumstcl- 
lar dust must be regenerated, presumably through grind- 
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ing collisions of planetesimals. Surveys for excess emis- 
sion in the infrared have demonstrated that at least 15% 
of nearby main sequence stars host debris disks (Habing 
et al. 2001; Bryden et al. 2006). As the sensitivity of in- 
struments improves, that fraction is steadily increasing. 
Already the Herschel astronomical observatory has iden- 
tified several new debris disk systems (e.g., Eiroa et al. 
2010, 2011; Thompson 2010), with more likely to be an- 
nounced soon as key programs wind down. Spatially re- 
solving the emission from the extended, faint dust disks 
is challenging, however, and so far only a handful of sys- 
tems out of hundreds of candidates have been spatially 
resolved. 

Resolved observations of debris disks reveal a wide 
range of morphologies, including broad disks, narrow 
rings, eccentricities, warps, and brightness asymmetries 
(e.g. Holland et al. 1998; Schneider et al. 1999; Wyatt 
et al. 1999; Heap et al. 2000; Kalas et al. 2005). Several 
of these morphological features have now been linked to 
the presence of planets in the disk; in fact, two of the 
three (so far) directly imaged extrasolar planetary sys- 
tems were predicted based on the properties of their dust 
disks (Fomalhaut and /3 Pictoris; Kalas et al. 2008; La- 
grange et al. 2010). That said, the clumpy structure 
predicted to result from large dust grains tracing plane- 
tary resonances has a more checkered history. Theoret- 
ical predictions for the magnitude of dust concentration 
in orbital resonances vary widely, with some groups pre- 
dicting a pronounced contrast in the surface density of 
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large particles (e.g., Ozernoy et al. 2000; Wyatt 2003, 
hereafter W03) and others predicting very little struc- 
ture even in the presence of massive planets (e.g., Kuch- 
ner & Stark 2010). A complicating factor is that disk 
structure is predicted - and, indeed, observed - to vary 
strongly as a function of wavelength. For example, sev- 
eral debris disks have now been observed to host smooth, 
extended haloes of small dust grains at radii far larger 
than the outer disk edges seen at longer wavelengths (e.g. 
Su et al. 2005, 2009). 

Millimeter-wavelength observations are preferred for 
tracing dynamical interactions with planets, since they 
trace the largest accessible dust grain populations that 
are the least sensitive to the smoothing effects of stel- 
lar radiation (Wyatt 2006). Early observations were 
promising. The Vega debris disk, in particular, exhib- 
ited a possible azimuthal asymmetry in the first 860 /im 
James Clerk Maxwell Telescope (JCMT) map of the 
disk structure (Holland et al. 1998), which was sup- 
ported by interferometric observations with the Owens 
Valley Radio Observatory (OVRO; Koerner et al. 2001) 
and the Plateau de Bure Interferometer (PdBI; Wilner 
et al. 2002). This result was followed by other claims 
of clumpy structure at millimeter wavelengths, includ- 
ing a suggestive observation of e Eridani in which the 
clumps appeared to move at super-Keplerian velocities 
between epochs, as expected for a resonant orbital con- 
figuration (Greaves et al. 2005). However, it was later 
pointed out that background galaxies may have con- 
tributed to the apparent asymmetry and rotation rate 
in e Eridani (Poulton et al. 2006), and observations at a 
wavelength of 350 /im did not confirm the presence of 
the 850 /^m clumps (Backman et al. 2009). More re- 
cently, clumpy structure has been claimed (and plane- 
tary properties inferred) in millimeter-wavelength maps 
of the disks around the young solar analogue HD 107146 
(Corder et al. 2009) and the multiple-planet host star HR 
8799 (Patience ct al. 2011). Unfortunately, more sensi- 
tive follow-up observations did not recover the claimed 
clumpy structure (Hughes et al. 2011). At the time 
of writing of this article, the observational evidence for 
clumpy millimeter-wavelength structure in debris disks 
is shaky at best. However, it should be noted that 
all millimeter-wavelength observations of debris disks to 
date have been plagued by low signal-to-noise, due to 
the tenuous nature of the dust emission (which rarely 
amounts to more than an earth mass of material). It is 
likely that sensitive future observations, e.g. with the 
Atacama Large Millimeter Array (ALMA), will uncover 
any resonant structure at lower contrast than the all-or- 
nothing levels that are currently observable. 

The Vega debris disk provides an interesting case - and 
is in many ways the prototype - of the study of clumpy 
structure in debris disks. The star itself is, of course, re- 
markably prominent both in the astronomical literature 
and in the night sky; it is an AOV star viewed nearly 
face-on to its rotation axis (Gulliver et al. 1994), located 
at a distance of 7.7 pc from the Sun (Perryman et al. 
1997; van Leeuwen 2007). Its age is estimated at 350 Myr 
(Song et al. 2000), and it is host to the first extrasolar de- 
bris disk ever observed (Aumann et al. 1984). Attempts 
to infer the disk size from IRAS observations resulted 
in estimates of --20" (150 AU; Harvey et al. 1984; van 
der Bliek et al. 1994); subsequent ISO data suggested 



that the disk may be even larger, up to 36" (280 AU) 
at a wavelength of 90 /xm (Heinrichsen et al. 1998). The 
initial reconnaissance of disk morphology took place in 
the millimeter, with apparent clumpy structure on 9" 
(70 AU) scales attributed to the presence of an orbit- 
ing Neptune-mass planet (Holland ct al. 1998; Koerner 
et al. 2001; Wilner et al. 2002, W03). Near-infrared 
interferometry with the Palomar Testbed Interferome- 
ter revealed a dust disk component much closer to the 
star, at a distance of only 4 AU, recently confirmed with 
IOTA (Ciardi et al. 2001; Defrere et al.'2011). After the 
flurry of work on the millimeter structure, it was with 
some surprise that Su et al. (2005) reported azimuthally 
symmetric dust emission observed with the Spitzer space 
telescope extending out to hundreds of AU at mid-IR 
wavelengths. The smooth structure was confirmed by 
subsequent Herschel observations at wavelengths of 70- 
500 /im; while the ring structure was resolved only at the 
shortest wavelengths, once again no azimuthal asymme- 
try was evident (Sibthorpe et al. 2010). Meanwhile, Gal- 
tech Submillimcter Observatory (CSO) observations at 
wavelengths of 350 and 450 /im resolved the ring struc- 
ture and reported tentative evidence of clumpy structure 
in the disk (Marsh et al. 2006). However, the position of 
the clumps in the CSO maps did not correspond with the 
clumps reported in the interferometric data; the authors 
suggested that perhaps the shorter- wavelength data were 
tracing 4:3 or 3:2 resonances rather than the 2:1 reso- 
nances detected at 860 /im and 1.3 mm wavelengths, and 
that two different grain populations may be responsible 
for the emission. However, the interferometric images 
were recently called into question by Pietu et al. (2011), 
who re-observed Vega with PdBI in a mosaicked obser- 
vation at a factor of two higher sensitivity than before 
and did not recover any clumpy structure. 

We present three new sets of observations of the Vega 
system at millimeter wavelengths, with the intention of 
helping to resolve the puzzling and contradictory obser- 
vations and models. The data were obtained with the 
Submillimcter Array (SMA)^^, the Combined Array for 
Research in Millimeter- wave Astronomy (CARMA)^^, 
and the Green Bank Telescope (GBT)^**, at wavelengths 
of 880 /im, 1.3 mm, and 3.3 mm, respectively. We de- 
scribe the observations in Section 2, present the results 
in Section 3, and analyze the data in the context of sev- 
eral historically relevant classes of models in Section 4. 
We discuss the results and summarize their significance 
in Section 5. 

2. OBSERVATIONS 

^■^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

^^ Support for CARMA construction was derived from the Gor- 
don and Betty Moore Foundation, the Kenneth T. and Eileen L. 
Norris Foundation, the James S. McDonnell Foundation, the As- 
sociates of the California Institute of Technology, the University 
of Chicago, the states of California, Illinois, and Maryland, and 
the National Science Foundation. Ongoing CARMA development 
and operations are supported by the National Science Foundation 
under a cooperative agreement, and by the CARMA partner uni- 
versities. 

^■^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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TABLE 1 
Basic Observing Parameters 



Date 



225 GHz opacity Track length (hr) RMS Sensitivity (in Jy beam ) 



341 GHz SMA - 4 GHz bandwidth 



2007 Apr 20 
2007 Apr 21 
2007 Apr 22 
2007 Apr 23 



0.05 

0.05-0.06 

0.04 

0.03 



2.6 
2.6 
2.2 
2.0 



227 GHz CARMA - 8 GHz bandwidth 



2011 Jan 20 


0.15 


3 


2011 Jan 23 


0.10 


5 


2011 Jan 30 


0.25 


4 


2011 Feb 2 


0.07 


2 


2011 Feb 7 


0.15 


3 


2011 Feb 8 


0.14 


6 


2011 Feb 9 


0.13 


6 



0.54 
0.53 
0.87 
0.87 
0.71 
0.61 
0.51 



''Angular resolution using natural weighting of visibilities. 

Data were collected in a series of 8 observing sessions 
on the GBT in spring 2009 and spring 2010. The GBT's 
90 GHz bolometer array, MUSTANG (Dicker et al. 2009) 
was used for these observations, giving a 9" FWHM beam 
on the sky. A total of 9.4 hours of useful integration 
time were acquired on the map of Vega. A 1.5x1.5 ar- 
cmin map around Vega was obtained with alternating 
Lissajous patterns, some centered at the position of the 
star and four others offset by ±15" in both north-south 
and east-west directions. Primary flux calibration was 
with reference to the asteroid Ceres- 1, using a light curve 
and brightness temperature provided by T. Mueller (pri- 
vate comm). Every 30-45 minutes the GBT beam, focus, 
and pointing was checked by a short (1 minute) obser- 
vation of a compact pointing calibration source within 
15 to 20 degrees of Vega. Pointing corrections were ap- 
plied offline in the data analysis. If the beam or focus 
showed evidence of drift from their optimal values, they 
were corrected by collecting a 15-minute Out-of- focus 
(OOF) holography map, analyzing the data on the fly, 
and applying the corrections to the telescope (Nikolic 
et al. 2007). Noise levels were established by means of 
a suite of Monte-Carlo bootstrap simulations in which 
the bolometer data were shifted by random amounts in 
time with respect to the antenna position data. Images 
were made by a series of IDL routines developed to ana- 
lyze MUSTANG data. These routines remove a common 
mode from all detectors, which is highly effective at elim- 
inating large-scale (>40") signals due to fluctuations in 
the atmosphere emission and cryogenic instabilities. A 
low-order polynomial is fitted to each detector to remove 
residual drifts; the individual detector weights are cal- 
culated from the RMS of the residual timestream; and 
the resulting data griddcd onto the sky. More details 
about the data reduction and calibration procedures can 
be found in Mason et al. (2010) and Dicker et al. (2009). 

The SMA data were collected during four eight-hour 
tracks on 20-23 April 2007 in the subcompact config- 
uration, which provided baseline lengths of 10 to 45 m 
between the eight six-meter antennas (see Table 1 for 
basic information about individual tracks) . The weather 
was excellent, particularly on the night of April 23. The 
225 GHz opacity was stable between 0.04 and 0.06, with 
the exception of a few spikes in humidity on the night 
of April 21. The correlator was configured for maximum 
continuum sensitivity, with uniform spacing of 32 chan- 



nels in each 104 MHz correlator chunk across the 2 GHz 
bandwidth in each sideband. The LO frequency was set 
to 340.8 GHz, placing the CO(3-2) line in the upper side- 
band. Flux calibration was carried out using Titan and 
MWC 349. Because of the range in right ascension cov- 
ered by the target and calibrators, atmospheric and in- 
strumental gain calibration was carried out using both 
MWC 349 and 3c345, with the former dominating the 
early part of the track and the latter the late part of the 
track. The derived flux of MWC 349 was 2.45 Jy; 3c345 
varied between 2.29 and 2.51 Jy over the four nights. The 
quasar J1801-I-440 was included in the observing loop to 
test the quality of the phase transfer to Vega, although 
the photosphere of Vega was also detectable on all four 
nights and served as a secondary check on the calibration. 
Errors in the nominal baseline solutions were evident in 
the quasar phases, which did not track each other well. 
However, we were able to derive corrected baseline solu- 
tions using only the calibrators present in the data set. 
The updated baseline solutions resulted in a 10% im- 
provement in signal-to-noise on the detection of Vega's 
photosphere. Routine calibration tasks were carried out 
using the MIR^^ software package, while imaging and de- 
convolution were accomplished using the MIRIAD soft- 
ware package. A naturally weighted image of the four 
data sets yields an rms noise of 0.94mJybeam^^ in a 
5'.'3x4'.'7 beam at a position angle of -80°. 

The CARMA data were collected in a series of seven 
short (3-6 hour) tracks between 20 January and 9 Febru- 
ary 2011 (see Table 1). The array was in its most com- 
pact "E" configuration, with projected baseline lengths 
between 5 and 65 m. The shortest baselines were pro- 
vided by the nine 6 m diameter antennas, while the 
longest baselines were preferentially between the six 10 m 
diameter antennas. The correlator was configured for 
maximum bandwidth and therefore continuum sensitiv- 
ity, with eight 500 MHz bands in each sideband for a total 
bandwidth of 8 GHz. The LO frequency was 227.5 GHz 
and the IFs of the bands were distributed between 2.3 
and 7.3 GHz in such a way as to avoid atmospheric ozone 
lines. The observing pattern was a three-point mosaic 
centered on the star position and the positions of the 
two clumps reported in Wilner et al. (2002). Obser- 
vations of Vega alternated with the quasar J1848+323, 
which was used to calibrate the atmospheric and in- 

^^ See http://cfa-www.harvard.edu/~cqi/mircook.htnil 
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Fig. 1. — CARMA maps of the 227 GHz continuum emission 
from the Vega system. The star symbol marks the position of 
Vega; thermal emission from the stellar photosphere is detected at 
the 6o" level. The top panel shows an image with naturally weighted 
visibilities, resulting in a beam size of 4'.' 2x2'.' 5 with position angle 
-88° . The bottom panel has been tapered in the visibility domain 
to match the spatial resolution of Figure Ic of Wilner et al. (2002), 
resulting in a beam size of 5'.'4x4'.'7 at a position angle of 32°. 
The contours are [2,3,4,5]x0.28mjy beam^ (the rms noise). The 
asterisk symbols (along with the star symbol) mark the pointing 
centers of the 3-point mosaic, and are chosen to coincide with the 
position of the emission peaks reported in Wilner et al. (2002). No 
thermal emission from the dusty disk is detected, despite shorter 
baselines and lower rms noise than the PdBI data set. 

strumental phase variations. No test quasar was in- 
cluded in the observing loop since Vega's photosphere 
was bright enough to be detected in most of the tracks 
(all but the two with the poorest observing conditions) 
and provided an excellent test of the quality of the 
phase transfer. The quasar 3c454.3 was used as the 
passband calibrator, and Neptune as the flux calibrator. 
The derived flux of J1848-h323 varied between 0.45 and 
0.63 Jy, yielding a photospheric flux for Vega of approx- 
imately 0.7mJy beam^^. Because Vega was a daytime 
source in January and the pointing on the 6 m anten- 
nas is poorer during the day, optical pointing was car- 
ried out every two hours, using the photosphere of Vega 



as the pointing calibrator (Corder et al. 2010). Cali- 
bration, imaging, and deconvolution were accomplished 
using the MIRIAD software package. The rms noise in 
the naturally weighted image combining all the data was 
0.28mJybeam~^ in a 4'.'2x2'.'5 beam at a position angle 
of -88°. 

3. RESULTS 

We detect the stellar photosphere in each of the three 
data sets, but only in the GBT data is there a marginal 
(3cr) detection of thermal emission from the dust disk. 
The disk emission is not detected in either of the two 
interferometric data sets despite lower rms noise and 
shorter baselines than previous observations that appar- 
ently showed clumpy dust emission. 

Fig. 1 shows the results of the CARMA observations. 
The stellar photsphere is independently detected in most 
of the seven short tracks, and the photospheric detec- 
tion for the combined data set is at the 6a level. The 
top panel of Fig. 1 shows the combined data set imaged 
with natural weighting, while the bottom panel shows the 
same visibilities imaged with a Gaussian taper to recover 
the lower spatial resolution with which the clumps were 
previously claimed to be seen. The asterisk symbols in 
the two panels mark the pointing centers of the 3-point 
mosaic (with the third point located at the star posi- 
tion) , chosen to coincide with the location of the putative 
clumps in the Wilner et al. (2002) data. No emission is 
detected at those locations in the CARMA data set de- 
spite a factor of two lower rms noise. The observations 
also have lower rms noise and richer (m,u) coverage than 
the OVRO data obtained by Koerner et al. (2001) using 
the 10.4 m antennas that have since been incorporated 
into CARMA. 

The SMA observations in Fig. 2 tell a similar story. 
The photosphere is pointlike, implying successful cali- 
bration of the instrumental and atmospheric gains, and 
is strongly detected at the 7a level. The higher fre- 
quency of the SMA observations requires flux scaling for 
comparison to the 1.3 mm data. Assuming a millimeter- 
wavelength spectral index of 2.8 for the disk, consis- 
tent with the long- wavelength Herschel Quxes reported in 
Sibthorpe et al. (2010), the expected integrated 880 /im 
flux of the emission peaks reported in (Wilner et al. 2002) 
is 20 and 12mJy for the northeast and southwest peak, 
respectively, with peak fluxes at the resolution of the 
data of approximately 7.4 mJy. Given the 0.9 mJy rms 
noise in the data, we would expect to detect emission 
from the clumpy dust disk at the 8a level (similar to the 
detection level of the stellar photosphere). However, no 
such emission was observed. This is comparable to the 
limits set by the CARMA data (Fig. 1) and the recent 
updated observations from the Plateau de Bure Interfer- 
ometer (Pietu et al. 2011). 

Taken together, these three interferometric data sets, 
observed with three different instruments with broad 
{u,v) coverage and uniformly low rms noise, place strong 
constraints on the amount of concentrated, clumpy emis- 
sion that might be present in the Vega debris disk. They 
also demonstrate that previous detections of clumpy 
emission from the Vega system using interferometers 
must have been spurious. The noise may have been diffi- 
cult to characterize near the edges of the primary beam, 
a problem mitigated in the recent data sets by the use of 
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Fig. 2. — SMA map of the 341 GHz continuum emission from 
the Vega system. Thermal emission from the stellar photosphere 
is detected at the la level. The visibilities have been imaged using 
natural weighting, resulting in a beam size of 4'.'7x5'.'3 with posi- 
tion angle -80°. The contours are [2,3,4,...] x0.94mjybeam~^ (the 
rms noise). Symbols are as in Fig. 1. No thermal emission from 
the dusty disk is detected, consistent with the CARMA (Fig. 1) 
and PdBI results (Pietu et al. 2011). 




Fig. 3.— 90 GHz GBT/MUSTANG map of continuum emission 
from the Vega system. The GBT map has been smoothed to 10" 
resolution (with 2" pixels) , indicated by the beam in the lower left 
corner. The rms noise varies with position due to the scan pattern, 
with the lowest rms closest to the pointing center, so the data are 
presented as a signal-to-noise map to emphasize the significance of 
the emission detection as a function of position on the sky. The 
emission is centered on the star position (peak signal-to-noise ratio 
of 4.7), and the integrated flux is 690itl30/xjy within a 16" radius 
from the star, indicating a marginal (3(t) detection of emission 
above the 310 /.tjy photosphere. 

mosaicked imaging (CARMA, PdBI) and smaller anten- 
nas (SMA). 

Given the smoothness of the disk emission indicated 
by the new interferometric observations, single-dish data 
are more desirable for characterizing the amount and lo- 
cation of emission from the dusty disk. The 90 GHz GBT 
map is presented in Fig. 3. The star position accounts 
for proper motion and is offset from the map pointing 
center by a few arcseconds. The location of the emis- 
sion peak coincides with the star position, although it is 



only determined to within 3" or so due to pointing uncer- 
tainties. The map has been smoothed to 10" resolution 
(with 2" pixels) . Because of the scanning pattern used to 
minimize atmospheric effects, the rms noise varies with 
position; it is lowest near the map center and highest at 
large distances from the map center. Fig. 3 shows a ±20" 
box around the star position. The emission peaks near 
the star position with a maximum signal-to-noise ratio 
of 4.7. The integrated flux within a 16" radius of the 
star is 690±130/iJy. To determine the expected pho- 
tospheric contribution at this wavelenghth, we use the 
most recent Kurucz model photosphere of Vega (R. Ku- 
rucz, priv. comm.), which is scaled to match the flux of 
7.17 Jy at 23.68 /im derived by Rieke et al. (2008). We ex- 
trapolate to the millimeter regime assuming that the flux 
scales as v^ in the Rayleigh- Jeans tail of the distribution. 
This yields photospheric fluxes of 4.4 mJy, 2.1mJy, and 
310 ^Jy at wavelengths of 880 /im, 1.3 mm, and 3.3 mm, 
respectively. The GBT data therefore imply a marginal 
(3cr) detection of 375 //Jy of emission above the photo- 
sphere, presumably from the dusty disk around Vega. 
Assuming a characteristic dust temperature of 70 K (Su 
et al. 2005, 2006) and a standard 90 GHz dust opacity of 
0.9cm^g^^g^ (Beckwith et al. 1990), this corresponds to 
a mass of 2.3xlO~'^M0, in good agreement with more 
detailed SED-based models (Su et al. 2005). Figure 4 
plots the measured photospheric fluxes at 880/im and 

1.3 mm along with the Kurucz model photosphere ex- 
trapolation to the millimeter spectral regime. A 3.3 mm 
upper limit (3cr) on the photospheric flux from Wilncr 
et al. (2002) is also plotted. The gray shaded region of 
the plot represents an extrapolation of the total 860 \xTa. 
flux from the disk reported by Holland et al. (1998), as- 
suming a spectral index of 2.8±0.1 estimated from fit- 
ting the long-wavelength Herschel data (Sibthorpe et al. 
2010). The expected steUar contribution at 860 /xm is 

4.4 mJy, while extrapolation of the Herschel data predict 
a larger 860 //m flux by about the same amount (which 
is within the measurement uncertainty). We therefore 
assume that the disk flux at 860 /im is 46 mJy. The GBT 
3.3 mm flux is only marginally consistent with both the 
Kurucz model photosphere extrapolation and the PdBI 
upper limit on the photospheric flux at 3.3 mm (Wilner 
et al. 2002), supporting the conclusion that the disk con- 
tributes to the measured 3.3 mm flux. A substantially 
flatter millimeter spectral index would weaken the in- 
ferred detection of 3 mm emission from the Vega debris 
disk. 

Given the low spatial resolution and signal-to-noise 
in the GBT data, it is difficult to draw strong conclu- 
sions about the spatial distribution of the emisson, al- 
though the implied integrated 90 GHz emission from the 
dust disk within 16" (125 AU) is approximately 375 /tJy. 
This is consistent with previous observations from the 
JCMT by Holland et al. (1998), who report a peak flux 
of 17.3±3.0mJy in a 14" beam at a wavelength of 860 /im. 
Scaling the 860 /im flux to the 3.3 mm wavelength of 
the GBT observations predicts 380 /iJy of dust emis- 
sion within 14" of the star, using the spectral index of 
2.8±0.1 derived from the long- wavelength Herschel data 
(Sibthorpe et al. 2010). The integrated emission in the 
JCMT map was larger, 46±5mJy, but the increasing rms 
noise with distance from the GBT map center makes it 
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Fig. 4. — Millimeter-wavelength SED of Vega. Flux values from 
this work are plotted with triangle symbols. Error bars mark 1-cr 
intervals, including both the noise in the image and a 20% system- 
atic flux uncertainty. Error bars are not included on the 1.3 mm 
data to avoid excessive cluttering around the four closely spaced 
points. The telescope name, angular resolution, and year of ob- 
servation is indicated in the legend. Data are drawn from Holland 
et al. (1998), Kocrner et al. (2001), Wilner et al. (2002), and Pictu 
et al. (2011). The disk is resolved out in the interferometric data 
(all points except the 860 /im JCMT flux from Holland et al. 1998 
and the 3.3 mm GBT flux), which are sensitive only to the pho- 
tosphere. A 3.3 mm upper limit (3o") on the photosphere from 
Wilner et al. (2002) is also included. The gray area represents an 
extrapolation of the 46±5 mjy integrated 860 fira flux measured 
with the JCMT (Holland et al. 1998) to longer wavelengths, as- 
suming a spectral index of 2.8±0.1 estimated from fitting the long- 
wavelength Herschel data (Sibthorpe et al. 2010). An extrapolation 
of the 2005 Kurucz model photosphere (solid line) is also plotted. 

difficult to detect very extended emission, so we limit 
ourselves to considering emission only within ~125AU 
from the star. Also, given the position uncertainty in the 
GBT map, it is difficult to determine whether the emis- 
sion peaks at the star position. The 3" uncertainty makes 
the GBT map at best marginally consistent with the 9" 
offset of the peak emission in the JCMT map; however, 
the photosphere represents a much larger contribution to 
the total flux at 3.3 mm than at 860 /im, which is likely to 
shift the location of the emission peak closer to the cen- 
ter as the disk emission is convolved with the photsphere. 
To summarize, while the GBT data are consistent with 
the JCMT map, on their own they can neither confirm 
nor disprove the presence of non-axisymmetric emission 
in the Vega debris disk. 

4. ANALYSIS 

The data presented in this paper, along with recent 
PdBI observations from Pietu et al. (2011), represent a 
marked departure from previous observations that indi- 
cated clumpy emission from the Vega debris disk. We 
therefore seek to characterize the parameter space con- 
sistent with the new observations. We consider three 
categories of models: (1) the dynamical model of a reso- 
nant planetesimal population from W03, similar to that 
presented in Wilner et al. (2002) and designed to re- 
produce the JCMT map; (2) a smooth, axisymmetric 
disk model, motivated by ring-like structure observed at 
shorter wavelengths; and (3) a "toy" model of clumpy 
debris disk emission that can be distilled into a two- 
parameter space to constrain exactly how clumpy (or 



axisymmetric) the disk emission may be while maintain- 
ing consistency with the data^^. We compare the data 
with these models in sequence in the following sections. 
While the GBT data are not sensitive enough to strongly 
constrain the degree of axisymmetry of the disk emis- 
sion at large distances from the star, they are useful for 
constraining the radial size of the flux distribution. We 
therefore include these data in the subsequent analysis. 

4.1. Resonant Clump Model 

Since the W03 model was generated to reproduce the 
clumpy emission reported in the older data sets (Holland 
et al. 1998; Koerner et al. 2001; Wilner et al. 2002), it 
is likely that it will be inconsistent with the more sensi- 
tive upper limits on clumpy disk emission from the PdBI 
(Pietu et al. 2011), SMA, and CARMA (this work). Nev- 
ertheless, for completeness we present a comparison of 
the new SMA, CARMA and GBT data with the W03 
model, scaled to the appropriate wavelengths. 

Figure 5 shows a comparison of the W03 model and 
(noiseless) simulations of each of the three data sets pre- 
sented in this paper. To simulate the interferometric 
observations, we sample the model at the same spatial 
frequencies as the data, then invert the visibilities and 
deconvolve the image with the same parameters as were 
used for the data, including the flux cutoff in the CLEAN 
procedure. The negative contours in the interferometric 
images are artifacts of the CLEAN algorithm applied to 
the data. The W03 model of the underlying planetesi- 
mal distribution is scaled to an appropriate flux at each 
wavelength. At 341 GHz, the photosphere is assumed to 
be 4.4 mJy and the total flux in the disk model plus pho- 
tosphere is scaled to match the integrated flux of 46 mJy 
reported for the JCMT detection (HoUand et al. 1998). 
At 227 GHz, the 850 /zm model is scaled in flux assum- 
ing a spectral index of 2 for the star and 2.8 for the 
disk emission. The 90 GHz model prediction is gener- 
ated by assuming a 310 /iJy photosphere and scaling the 
flux within 16" of the star position to the integrated flux 
in the GBT map. The resulting simulated map demon- 
strates that the GBT data are not sensitive enough to 
rule out a resonant emission model. While the SMA 
data are similarly consistent with the W03 simulations, 
the CARMA data are clearly inconsistent with the model 
prediction. 

4.2. Axisymmetric Model 

Since a clumpy model is inconsistent with the observa- 
tions, it is desirable to characterize the parameter space 
allowed for a smooth, axisymmetric model of the disk 
emission. A smooth model is consistent with shorter- 
wavelength observations of the Vega debris disk (e.g., Su 
et al. 2005) and with submillimeter observations from 
Marsh et al. (2006), and is perhaps the most natural 
explanation for the lack of detectable millimeter-wave 
emission in the interferometric data sets. We therefore 
investigate the parameter space allowed by the combina- 
tion of SMA, CARMA, and GBT data presented in this 
paper. 

^^ We also tried an eccentric ring model motivated by the double- 
peaked structure in the JCMT and GBT maps, but it required 
emission too concentrated to be consistent with the SMA and 
CARMA data 
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Fig. 5. — Simulated observations of the W03 resonant model (left) with each of the different observing setups presented in this paper. 
If the model were correct, we would expect to detect dust emission at the 4a level with the SMA at 341 GHz (center left) and at the 
6(T level with CARMA at 227 GHz (center right). The 90 GHz GBT observations are consistent with the W03 model primarily because 
they do not have the spatial resolution or sensitivity necessary to discriminate between this and other models. The observations are 
simulated without noise and with identical imaging parameters and flux cutoffs to the images in Figs. 1 and 2 (the negative residuals in 
the interferometric maps are artifacts of the CLEAN algorithm applied to the data) , but the W03 model is clearly inconsistent with the 
observations, particularly the CARMA data. Symbols are as in Figs. 1-3, including SMA contour levels of [2,3,4,...]x0.94mjy beam~^ and 
CARMA contour levels of [2,3,4,...]x0.28mjy beam-^ 

We assume a smooth, axisymmetric disk with inner 
radius i?in and width AR (so that the disk outer ra- 
dius i?out = -Rin + AR). We assume that the disk is 
viewed face-on {i—Q). The inchnation of the stellar ro- 
tation axis to our line of sight is closer to 5° (Aufden- 
berg et al. 2006), but the true disk inclination is not 
well determined. Since many previous papers modeling 
the Vega debris disk have assumed a face-on orientation 
for simplicity, we follow suit; a difference of a few de- 
grees does not alter our conclusions. We assume constant 
surface density between Ri^ and i?out, with temperature 
decreasing as T{R) ex. R~^-^ and opacity constant with 
radius. Effectively, this amounts to an assumption that 
flux also decreases with distance from the star as R~^-^. 
We make no assumptions about dust grain composition, 
emission or absorption efficiency, or opacity. The tem- 
perature and mass of the dust are thoroughly degenerate, 
and we concern ourselves only with the distribution of 
flux in the disk, i.e., the product of surface density, tem- 
perature, and opacity. For standard assumptions about 
the dust grain composition and opacity, the tempera- 
tures assumed in our model are consistent with those in 
the hterature (e.g., Su et al. 2005). We further adopt 
the observational constraints that the total flux of star 
and disk at 860 /im is 45.7mJy (Holland et al. 1998); 
the photospheric flux at 341 GHz is 5.7 mJy; the flux at 
227 GHz scales from 860 /im with a spectral index of 2.8; 
and the total disk flux within 16" from the star at 90 GHz 
is 375 /xJy. Under these constraints, we vary i?in and AR 
to generate a grid of models to be compared with the 
data. For the interferometric data sets, we generate a 
model of the disk only (without the photosphere), and 
sample the visibilities with the same spatial frequencies 
as the observations using the MIRIAD task uvmodel. We 
then compare the synthetic visibilities with the observed 
visibilities, real and imaginary, and calculate a x^ value. 
For the GBT data, we create a model image of the photo- 
sphere and disk, convolve it with an appropriately sized 
beam, and calculate the reduced x^ value in a 16" square 
region around the star position, assuming the 75 /iJy rms 
noise in the map center derived from simulated observa- 
tions in Section 2. 

Fig. 6 presents the results of the simulations, and 
marks off the parameter space allowed by the SMA, 
CARMA, and GBT observations, given the aforemen- 



FlG. 6. — Confidence interval contour plot showing the regions of 
parameter space for a smooth disk model allowed by the combina- 
tion of SMA (5" resolution), CARMA (4"), and GBT (10") data. 
The abscissa corresponds to the inner radius of an axisymmetric 
disk while the ordinate corresponds to the width of the disk (so 
that the disk outer radius -Rout = Rin + A_R). Small inner radii 
and small disk widths are ruled out by the interferometric data, 
while large inner radii are ruled out by the GBT data. 

tioned assumptions. Small values of i?in and AR are 
ruled out by the interferometric data, since emission con- 
centrated too close to the star position or in too nar- 
row a ring around the star should have been detected 
in the CARMA or SMA data sets. Models with larger 
Rin and AR exhibit smoother emission that would be 
undetectable by the interferometric observations; how- 
ever, based on the GBT detection of 375 /im within 16" 
(125 AU) of the star, Rin must be smaller than this aper- 
ture. The flux distribution in the GBT data places con- 
straints on smaller radii as well. In the context of a 
smooth, axisymmetric disk model, an inner radius of 
-20-100 AU is preferred, with a width >50AU. Inner 
radii <15 AU and >110 AU are ruled out by the CARMA 
and GBT data, respectively, while widths <30AU are 
ruled out by the interferometric data. 

4.3. Clumpy Toy Model 

Although the data presented in Pictu et al. (2011) and 
this work arc consistent with an azimuthally symmetric 
distribution of emission, it is still true that some mil- 
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limeter data sets (Holland et al. 1998; Marsh et al. 2006) 
exhibit marginally significant evidence for a clumpy dis- 
tribution of emission. In an attempt to quantify exactly 
how clumpy the Vega debris disk may be given the avail- 
able data, we consider a toy model consisting of a smooth 
component underlying two Gaussian "clumps" of vari- 
able FWHM. We make several simplifying assumptions: 
(1) the location and relative brightness of the two clumps 
are fixed at 9" northeast and southwest of the star posi- 
tion, respectively, with the southwest peak 0.3 times that 
of the northwest peak, consistent with a literal interpre- 
tation of the JCMT data; (2) the inner and outer radius 
of the smooth component are fixed at 60 and 200 AU, re- 
spectively, to maintain consistency with the results of the 
axisymmetric modeling from Section 4.2; and (3) the sur- 
face density of the smooth component is held constant, 
while the temperature is assumed to decrease with dis- 
tance from the star R as R~^'^. As in Section 4.2, we 
also normalize the 880 /xm flux to the measured JCMT 
value from Holland et al. (1998), assume a spectral in- 
dex of 2.8 to extrapolate to 227 GHz, and normalize the 
90 GHz disk flux to the measured 375 /xJy within 125 AU 
from the star. Given these assumptions, the only vari- 
ables are the fraction /dump of total flux in the clumpy 
component (with the rest distributed across the smooth 
component), and the FWHM of the Gaussian clumps. 
While we do fix the position of the clumps in the disk 
to match the JCMT observation, the result will be more 
generally applicable since the interferometric data would 
be capable of detecting comparably clumpy structure at 
virtually any position in the map. 

We first examine only the constraints from the inter- 
ferometric data, since these data are the most sensitive 
to the presence or absence of flux concentrations. Fig. 7 
shows the parameter space permitted by the SMA and 
CARMA observations: as expected, a smooth bright- 
ness distribution is favored over a clumpy distribution, 
although very diffuse and extended clumps (with FWHM 
> 100 AU) are permitted. It should be noted that the 
largest size scale we consider is ^120 AU, since this is 
roughly the largest size scale to which the data are sensi- 
tive, given the ~1 mm wavelength of observation and the 
~10m baselines in the CARMA and SMA data. How- 
ever, even when diffuse clumps are included, they should 
not account for more than about half of the total flux 
in the image to maintain consistency with the interfer- 
ometric data. The constraints from the GET data are 
somewhat orthogonal: due to the flux concentration in 
the southwest corner of the GET map, the data favor 
models with a clumpy underlying flux distribution. The 
overlapping preferred parameter space for all three data 
sets occurs near the upper right corner of the plot, with 
large, diffuse clumps accounting for approximately half 
of the flux in the data, and the rest distributed in an ex- 
tended smooth component. However, considering the low 
signal-to-noise ratio for the extended (non-photospheric) 
flux in the GET map, we note that the constraints from 
the interferometric map are the strongest, and require 
the smooth component of the flux distribution to domi- 
nate. The CARMA and SMA therefore demonstrate that 
most (^50%) of the flux is distributed in a smooth, ax- 
isymmetric ring around the star, although a small degree 
of concentration on very large scales is permitted. This 
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Fig. 7. — Confidence interval contour plot showing the regions 
of parameter space for a clumpy disk model allowed by the inter- 
ferometric (SMA and CARMA) data. The abscissa corresponds to 
the FWHM of the Gaussian clumps in AU, while the ordinate cor- 
responds to the fraction of flux in the clumpy rather than smooth 
component. 

result does not depend strongly on the assumed values 
of Rin and AR. 

5. DISCUSSION 

Modeling of the three new data sets favors a smooth, 
axisymmetric flux distribution for the Vega debris disk, 
with an inner radius of 20-100 AU and a broad radial 
width of at least 50 AU. Diffuse emission on large scales 
may account for up to about half of the flux, but the ma- 
jority must be smoothly distributed. We have considered 
only a small subset of possible flux conflgurations for the 
disk, but the results are more generally applicable. In 
particular, the clumpy toy model demonstrates that any 
flux concentration on spatial scales of less than '-^lOO AU 
cannot dominate the millimeter-wavelength flux distri- 
bution. While we chose a two-clump model to approxi- 
mate the morphology of the JCMT data, the result can 
be applied to any concentrated flux distribution. If the 
flux were concentrated into a single clump, for exam- 
ple, it would account for even less than half of the total 
flux. Nor is the location of the clumpy emission restric- 
tive. The interferometric data effectively rule out clumpy 
emission within the 30" primary beam of the SMA at 
880 /im and within the ~30" x45" approximately ellipti- 
cal region covered by the 3-point CARMA mosaic. 

It is unlikely that the clumpy structure detected by Ko- 
erner et al. (2001) and Wilner et al. (2002) could have 
redistributed itself in the ~ 10-year interval since those 
observations. Resonant patterns, like those described by 
W03, change and move on timescales comparable to the 
orbital period of the planet inducing the resonance. For 
the putative Neptune-mass planet located 65 AU from 
the star described in W03, the orbital period is approx- 
imately 300 yr. This is too long to account for a large- 
scale redistribution of material on a 10-year timescale 
(moreover, resonant structures are typically long-lived; 
they are not expected to dissolve on orbital timescales, 
much less fractions of an orbit). Stellar radiation forces 
can also alter the spatial distribution of dust; however, 
the grains responsible for the observed emission are likely 
millimeters in size, too large to be expelled by radiation 
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pressure, and the Poynting-Robertson drag timescale at 
the 70 AU radius of the claimed clumpy structure is of 
order 1 Gyr. The dust distribution could also change 
because of collisions: the millimeter-sized grain popula- 
tion could be removed by coUisional grinding, or be ren- 
dered invisible by growth into larger bodies with smaller 
millimeter- wavelength opacity. Both destruction and ag- 
glomeration occur on timescales at least as long as the 
coUisional timescale, which can be estimated as Porb/T, 
where Port is the orbital period and t is the vertical ge- 
ometric optical depth of the disk (e.g., Lagrange et al. 
2000). For a ring whose annular width is comparable to 
its radius, the vertical optical depth is roughly the ob- 
served ratio of infrared to bolometric luminosity, some- 
times denoted tjh (e.g., Backman 2004; Chiang et al. 
2009). Since the Vega debris disk exhibits Tiji—2xl0^'^ 
(e.g., Lagrange et al. 2000), the collision timescale must 
be orders of magnitude larger than the orbital timescale. 
Finally, one can imagine a scenario wherein the clumps 
observed ~10 years ago represent a transient burst of 
dust released from a catastrophic collision between larger 
parent bodies - dust whose surface brightness is today 
too small to be detected because of spreading of the 
ejecta. But even at extreme ejecta speeds of lkms~^, 
dust could only have spread across a few AU, which is 
an order of magnitude smaller than the linear scale cor- 
responding to the beam size of the observations; further- 
more, spreading by Keplerian differential rotation occurs 
over timescales at least as long as the local orbital pe- 
riod of 300 yr. In sum, it is difficult to imagine how 
the clumpy structure could have disappeared in the ~10 
years that have elapsed between observations. 

The interferometric data rule out millimeter- 
wavelength morphologies that concentrate the brightness 
on small scales; this includes a narrow belt or ring that 
would concentrate the flux radially. Ring widths of less 
than 30 AU are ruled out at the 3a level by the inter- 
ferometric data, and widths of >1G0AU are preferred. 
One emerging framework for understanding debris disk 
structure asserts that in some systems, large particles 
may be concentrated into a narrow "birth ring" in which 
collisions generate small grains carried onto broader 
orbits by effects of stellar radiation (Strubbe & Chiang 
2006). This theory has had notable success in explaining 
the differing radial extent of millimeter-wavelength and 
scattered light data, for example in the disk around /3 
Pictoris, which exhibits a ratio AR/R of <0.5 (Wilner 
et al. 2011). However, the Vega debris disk is one of sev- 
eral for which spatially resolved millimeter data require 
a broad band {AR/R > 1) rather than a narrow ring of 
flux (see also HD 107146 and HR 8799; Hughes et al. 
2011). At the coarse spatial resolution and limited sen- 
sitivity of current observations, not all systems fit neatly 
into a birth ring paradigm. Millimeter-wavelength maps 
display a variety of morphologies relative to scattered 
light. 

It should also be noted that the absence of significant 
clumpy structure on <100 AU spatial scales does not rule 
out the presence of planet-mass companions. There are 
several relevant points to consider. The W03 model only 
predicts the underlying planetesimal distribution, and 
does not include a detailed treatment of the effects of 
collisions on the spatial distribution of millimeter dust 
grains. There is some evidence that coUisional smoothing 



may wash out evidence of planetary resonances even for 
the large dust grain populations that dominate the emis- 
sion at millimeter wavelengths (Kuchner & Stark 2010), 
and that far more sensitive observations will be needed 
to observe the subtle contrast generated by dynamical 
interactions with orbiting planets. Furthermore, taken 
as a whole, the observational evidence so far points to a 
complex, multi-component debris disk around Vega. The 
data presented in this paper rule out a centrally concen- 
trated brightness distribution (which would indeed pro- 
duce a detectable flux concentration on <100 AU spatial 
scales) and instead favor an inner radius for the dust disk 
of 20-100 AU and a width >50 AU. However, the near-IR 
interferometry (Ciardi et al. 2001; Defrere et al. 2011) 
and a forthcoming re- analysis of the resolved images and 
SED taking into account all available mid- and far-IR 
data (K. Su, private communication) provide evidence 
for multiple dust belts in the system, including a warm 
asteroid belt at a radius of a few AU. The presence of 
multiple dust belts at different radii is suggestive of a 
planetary system. For example, the HR 8799 system is 
known to host a system of at least four planets, whose 
orbital radii are bracketed by warm and cold dust belts 
inferred from models of the system SED (Chen et al. 
2009; Suet al. 2009). 

The series of recent observations calling into question 
claims of clumpy structure in millimeter-wavelength de- 
bris disks (e.g., Pietu et al. 2011; Hughes et al. 2011, 
this work) also has unclear implications for future ob- 
servations. While there is little doubt that debris disk 
structures like warps, rings, and brightness asymmetries 
are tied to the presence of planets in the disk (e.g. Wy- 
att et al. 1999; Heap et al. 2000), it may be necessary to 
achieve substantially greater sensitivity levels than ini- 
tially thought to observe the density concentrations in- 
duced by orbital resonances in the large-grain dust dis- 
tribution. Fortunately, such observations are rapidly be- 
coming possible with the advent of the ALMA. The un- 
precedented sensitivity of ALMA combined with the re- 
cent discovery of planets orbiting in debris disks (Kalas 
et al. 2008; Marois et al. 2008; Lagrange et al. 2010) 
will provide an exciting testbed for determining the mag- 
nitude and wavelength dependence of density perturba- 
tions induced by planet-disk interactions in a range of 
different systems. The next-generation MUSTANG2 in- 
strument on the GBT also promises to be a valuable ad- 
dition to the suite of northern-hemisphere instruments. 
With 25 times greater sensitivity than the current MUS- 
TANG detector, the structure of the Vega debris disk 
will observable at high fidelity. MUSTANG2 will also be 
able to resolve a number of other nearby debris disks, 
opening up a new long- wavelength regime of debris disk 
imaging. 

6. SUMMARY 

New observations of the Vega system at a range of mil- 
limeter wavelengths support the observations of Pietu 
et al. (2011), indicating that previous detections of 
clumpy structure in the Vega debris disk were spuri- 
ous. We have analyzed new interferometric data sets 
from the SMA and CARMA, as well as a single-dish ob- 
servation with the GBT, in the context of several models 
of disk structure. We flnd that the interferometric data 
in particular require a disk with a large (20-100 AU) in- 



10 



Hughes et al. 



ner radius and a broad (>50 AU) radial width. Clumpy 
structure on <100AU scales accounts for less than half 
of the total 860 ^m flux observed by the JCMT (Holland 
et al. 1998). The disk must therefore be dominated by a 
smooth and largely axisymmetric millimeter flux distri- 
bution. 
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